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Abstract 


In this paper, we investigate the nature of the singularity in the spherically 
symmetrical, shear-free, gravitational collapse of a star with heat flux using a 
separable metric [||]. For any non-singular, regular, radial density profile for 


a star described by this metric, eq. (2.1), the singularity of the gravitational 
collapse is not naked locally. Our results here unequivocally support the 
Strong Cosmic Censorship Hypothesis. 
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I. INTRODUCTION 


The problem of whether the gravitational collapse leads, under general astrophysical 
circumstances, to a black hole or to a naked spacetime singularity - the Cosmic Censorship 
Hypothesis (CCH) - is one of the important problems of modern physics. However, this has 
proven to be a very difficult question to answer. It is therefore important for the proper 
mathematical formulation of this question that we study as many of the gravitational collapse 
situations as possible. It is expected that such studies would help sharpen the statement of 
the CCH. 

In a recent work []^ we have studied spherically symmetric, separable metric spacetimes 
with energy/heat flux. The source matter in those spacetimes can satisfy any equation of 
state, in particular, it is allowed to be of the form p = ap where p and p are the pressure and 
the density of the matter and a is a constant. We consider here the formation of singularities 
in the gravitational collapse of a star described using those spacetimes. 

The organization of this paper is as follows: In §|1| we outline the features of the metric 
of the shear-free, spherically symmetric, separable-metric spacetime with heat flux. This 
class of solutions is then analyzed in §|T| for the existence of the strong curvature spacetime 
singularity. In §|^ we argue that the solutions can model a star with heat flow and an 
extended atmosphere. In §[V] we analyze the radial light cone equation for this class of metrics 
and show that a horizon always forms before the singularity. We show that for a regular, 
non-singular, density profile of a star, the singularity is not even naked locally. Finally, in 
§[V|, we discuss the implications of our results for the Cosmic Censorship Hypothesis. 


II. THE SPACETIME METRIC 


The spacetime metric of all the spherically symmetric, shear-free, separable metric solu¬ 
tions with matter allowed to possess the equation of state p = ap is |]I| 


ds‘^ 


dt^ + 2 {y'Y dr'^ -|- {dO"^ -|- sin^ {6) d(j)‘^) 


( 2 . 1 ) 
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where y = y{r), R = R(t) and an overhead prime has been nsed to denote the derivative 
with respect to the radial coordinate r. The coordinates {t,r,6,(j)) are comoving. Note that 
the spatial part of the metric ( p.l|) is, in general, non-flat. 

The non-vanishing components of the Ricci tensor for the metric (p.l|) are 


Roo — ~3 -l- 


R R2 


( 2 , 2 ) 


Roi — 2 


Ry^ 

R y 


(2.3) 




(2.4) 


(| + 2 | 


( 2 . 6 ) 


Rrt-i = sin^ 9R' 


22 


( 2 . 6 ) 


The Ricci and the Kretschmann scalars of the metric ( |2.1|) are 

^ 6 fR R^\ 1 

7 ^ = — I - + 


1/2 \R Ry ym^ 


(2.7) 


K* _ ryabcd r) 

A_/ — it ■tvnhrd — 


3 - 4 R2 + 12 


R 12 


y^ R^ y^ R^ y^ R2 

Clearly, the metric and the spacetime are singnlar for y = 0 and/or R = 0. 

The flnid fonr-velocity f/“ and the flnid fonr-acceleration Ua = Ua-bU^ are given by 


( 2 . 8 ) 


= -(5“o 

y 


(2.9) 


f7,= (0,^,0,0) 

y 


( 2 . 10 ) 


where a semicolon denotes a covariant derivative. 
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We observe that the Rqi component of the Ricci tensor or, equivalently, the Gqi com¬ 
ponent of the Einstein tensor, is non-vanishing. This implies the presence of non-vanishing 
energy-flux in the spacetime under consideration. The energy-flux can be physically inter¬ 
preted in different possible ways. Firstly, there could be bulk motions of matter. Secondly, 
there could be “thermodynamical” heat flux arising from microphysical considerations of 
matter. For the purpose of a stellar application here, we shall consider the interpretation 
arising out of thermodynamical considerations of the stellar matter. 

As a result, the energy-momentum tensor of the matter fluid is of the form 

Tab = (p + p)UaUb + PQab + QaUb + QbUa (2.11) 

where p is the pressure, p is the density of the fluid and = (0, g, 0, 0) is the radial heat 
flux four-vector. Note that both shear and rotation vanish for the metric (p2.1| ). 

The expansion, 0, heat flux, g, density, p and pressure, p, are given by 

Q=U^-,a = -^ ( 2 . 12 ) 

y R 


-1 R 
y‘^y' R^ 


(2.13) 


P = 


R2 


+ 3 R^ 


(2.14) 




2R 


P = 


(2.15) 


2p2i?2 ^2^2 y2Jl 

Note that p oc \/y^ and that the radial function is not hxed by the held equations. Hence, 
we can, for a non-singular density distribution, choose a nowhere vanishing radial function 
y(r). Note, however, that for yir) = r or p' = 1, the initial density prohle is singular at 
r = 0 i. 

From eq. (|2.14D and (|2.15|), we obtain 


R = 


fR 

6 


2/2R2 


(P + 3p) 


(2.16) 
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It is then seen from eq. (p.l6|) that the equation of state for the collapsing matter closes 
the system of ordinary differential equations obtainable from the held equations and leads 
to dehnite dynamics for the collapsing matter. 


We note that, for the metric the equation of state for the collapsing matter could 

change with the progress of the collapse. In other words, we emphasize that the metric ( p.lD , 
although initially obtained for a specihc equation of state of the form p = ap for a being 
a constant, admits any equation of state for the collapsing matter. In fact, the equation of 
state is an additional input of a physical nature that is to be supplied in order to determine 
the dynamics of collapsing matter. 

Further, the change in the equation of state of the stellar matter could “halt” the collapse 
and may result in a semi-stable or stable matter conhguration. However, since our interest 
here is in the end-result of the unstoppable gravitational collapse, we, in what follows, 
assume that the gravitational collapse is unstoppable. Hence, our further results below 
apply to “unstoppable” gravitational collapse. 

We also emphasize that in the present description we could begin with an initial stellar 
density conhguration that is non-singular for all r. Of course, the singularity will develop as 
a result of the time evolution when the physical radius of the collapsing stellar shell vanishes 
at some moment of the time, ie. when i? = 0 at some t = to, say. Further, from symmetry 
considerations, such a singularity will, of course, be located at r = 0. 

III. THE SPACETIME SINGULARITY 

From the physical point of view, a singularity could be termed strong if a test body is 
crushed to zero volume as it approaches the singularity. Hence a sufficient condition for the 
singularity to be a strong curvature type as dehned by Tipler et al 0 is: 

For at least one null geodesic with an affine parameter s, with s = 0 at the 
singularity, the following holds 

lim s'^RabK^K^ > 0 (3.1) 


5 





where is the tangent vector to the geodesic in qnestion and Rab is the Ricci 
tensor. 


If this condition is satished then the spacetime is not extendible [Q. Snch spacetimes with 
strong curvatnre singnlarities are to be considered in view of the CCH. 

The geodesics for the metric (p^.l| ) are easily obtainable by considering the metric as 
a Lagrangian: 

2C = —yH'^ + + y‘^6^ + y‘^ sin^(6') 0^ (3.2) 

where an overhead tilde denotes a derivative with respect to the affine parameter s along 
the geodesic and 2C has valnes +1,-1 and 0 for the spacelike, timelike and nnll geodesics, 
respectively. 

The geodesic eqnations are then obtainable as 




dt 


(3.3) 


d 


- {RSy'f) = c. 


(3.4) 


A 

ds 


{R^y^O^ + cot^ 6 


= 0 


(3.5) 


d 


— (^R^y^ sin^ 6*0) = 0 


(3.6) 


where A; is a constant of integration from the 0-eqnation. 

Let ns then consider the radial nnll geodesics of the metric (|2.1|) for the pnrpose of the 
above condition since it reqnires the existence of at least one nnll geodesic for its pnrpose. 
Now, the fntnre-directed tangent to the radial null geodesic has the following components 


_ dt Ki 


K^ = — = 


ds y‘^R 


(3.7) 


K'- = -= 

ds 2yy'R‘^ 


(3.8) 
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where K 2 are constants. Since is null, we have K 2 ^ = 2Ki^. 
We now have 




(R 

+-h 2- , 

\R R^ 2y^R^ 


Kl 


+ 


2RK1K2 


[p + p) - -^y Ry 


ym 


y^ 


(3.9) 

(3.10) 


We can now, after some manipulations based on the energy conditions, in particular, {p + 
pY > a (see Appendix), show that 


lim s^RabK^K^ —^ 00 (3.11) 

Note that we have s = 0 at the singularity i? = 0. 

Hence, the singularity of the metric (p.l|) is a strong curvature singularity in the sense of 
Tipler et al. The singularity in a spherical gravitational collapse with heat flow considered 
here is then a strong curvature singularity always. 

This completes our description of the shear-free, spherically symmetric, separable metric 
solutions of the Einstein held equations that represent an inhomogeneous distribution of 
imperfect matter with heat hux. Different energy conditions can also be satished for the 
metric (13) • (See Appendix 0. This spacetime is also of interest to cosmology g].) 


IV. STELLAR DENSITY PROFILE 

Any initial density conhguration of the star must be non-singular for all r. For example, 
by noting that p oc l/y"^, the radial function |I| could be chosen as 

y = l+exp(^) (4.1) 

where Tq is the boundary radius and L is the thickness of the boundary layer. The corre¬ 
sponding star has density that decreases with r and at the boundary radius, Tq, there is a 
sharp fall in density within the layer L. Note that the density does not vanish outside the 
star but can be very low as compared to that inside the star. In the case of our sample 
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profile, eq. (|4.1|) , the star has an atmosphere of exponentially decreasing density. Essen¬ 


tially, the star is embedded in a cosmological spacetime of some sort and the spacetime is 
not necessarily asymptotically flat. The presence of an extended atmosphere enveloping the 
star has been observed with real stars, for example, with the Snn. We emphasize that we 
do not assume the profile 0 to obtain the results presented here. 

Also, the “exterior” of the star as considered here needs to be specified. If we consider 
the exterior of the star as being described by the Vaidya spacetime then we obtain the 
matching condition p as: 

pj: = V2yy'\j: (4.2) 


where S is the timelike hyper-snrface across which the interior and the exterior spacetimes 
are matched. This condition restricts the allowed forms of the temporal fnnction R{t) at the 
stellar bonndary. However, for onr pnrposes, it is not necessary to consider these restrictions. 

Moreover, note that q oc —R/y'. Dnring the collapse, we mnst have i? < 0 (from (|2.12|) ). 
Hence, with onr assnmption that y' > 0, the heat flux is positive, ie. heat flows from 
high temperature regions to low temperature regions. Our model of a star with heat flow 
therefore has a higher temperature at the center. The temperature decreases as we move 
away from the stellar center. This is what we expect in any realistic model of a star. In our 
model, the density will be a decreasing function of the radial coordinate r when y' > t] and 
this is what we shall assume. 

Of course, y{r) can be constrained by other considerations. For example, the thermal 
stability of the stellar structure, heat transport equation etc. will, in a manner similar to 
that of the Newtonian theory, yield some physical constraints on the radial distribution of 
density and, hence, on y{r). In fact, it is these considerations that must be used to obtain 
the detailed stellar model in a manner similar to that of the Newtonian theory. (However, 
we wish to point out here that these are not relevant to the subject of the present paper 
mainly because the present considerations apply to all the forms of the function y{r).) 


We should also keep in mind that we require the spacetime manifold to be locally flat 




in the neighborhood of the origin, r = 0, of the coordinate system, as it should be for any 
point of the spacetime manifold. In our case, a small circle of coordinate radius e with center 
at the origin has circumference of 27re. On the other hand, the circle has the proper radius 
y/2y'e. Then, requiring that the ratio of the circumference to the proper radius of the circle 
to be 27r in the neighborhood of the origin, we obtain the condition ~ l/v^- This 

condition must be imposed on y{r) to obtain a realistic density prohle of a star. 


V. OUTGOING RADIAL NULL GEODESICS 


The radial null cone equation is 


dr y 


( 6 . 1 ) 


where we have to choose the positive sign for the radially outgoing null geodesics. 

The physical radius of the collapsing shell of matter is Rph = yR- The collapsing, 
spherically symmetric shell of a star must have some physical radius Rph at f = to- The 
physical radius becomes zero at t = 0 signifying that the shell in question has collapsed 
to a singularity. The singularity will, of course, be at the radial coordinate r = 0 due to 
symmetry considerations. The singularity of the metric is then at f = 0 and r = 0. 

Further, we could choose y{r) = r, that is, y' = 1. However, as remarked earlier, this 
choice corresponds to a singular initial density prohle and, hence, is unacceptable to us 
since a normal star does not show such a density prohle. Therefore, we use the non-singular 
density prohle henceforth. 

To determine whether the singularity is globally visible or not, we need to look at the 
convergence of the principal null congruence and determine the condition for the formation 
of the outermost trapped surface. Now, the outgoing principal null vector for the metric 
( pTl) is 


1 

71 


+ -j }-— 

y \/2y'R 


5 “ 


(5.2) 
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To test for the existence of a trapped surface, we set the expansion rate for to zero and 
look for a real solution that is positive, Rah- This solution corresponds to the apparent 
horizon - the time history of the marginally trapped surface of the spacetime geometry of 


the metric (|2.1|) . 

Moreover, since we are dealing with separable metric functions, we expect to obtain a 
condition only on the time-dependent part of the metric, that is, R{t). The condition is 


R = 


1 

72 


(5.3) 


The trapped surface then always covers the singularity since R must attain the above value as 
it becomes (negatively) unbounded with the progress of the collapse. Hence, the spacetime 
singularity is not globally visible in all of the cases being considered below. 

When we consider a non-singular density prohle, we shall have 


lim - = L 

r—>0 y 

say, where to is the appropriate positive and hnite limiting value. Now, for the geodesic 
tangent to be uniquely dehned and to exist at the singular point, f = 0 and r = 0, the 
following must hold so that an outgoing, future-directed photon trajectory exists at the 
singularity: 


lim - = lim — = Xq (5. 

r ^ 0 r r ^ 0 dr 
t^O t^O 

where Xo is required to be real and positive. As we approach the singular point, we have, 

Xq = a/ 2 lim — Rit) = a/2 \im. to Rit) (5.5) 

r —> 0 y t^O 

Clearly, the limit in question implies that Xo = 0. Hence, there does not exist a real and 
positive tangent to the null geodesic at the singularity if the initial density prohle corresponds 
to a density distribution that is non-singular. 

Therefore, no null geodesics emerge from the singularity and communicate to an external 
observer. Then, for any non-singular density prohle of a star with heat hux the singularity 
of collapse is not visible to any observer. 
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VI. DISCUSSION 


A genuine spacetime singularity is a singular point of the very structure of the spacetime. 
The usual laws of classical physics break down at such a point of the spacetime manifold. 
It is therefore impossible to account for the effects of the singularity if it were to causally 
affect its exterior. Further, the usual laws of quantum physics, relativistic or non-relativistic, 
also fail to hold at such a spacetime location. This is mainly because such laws have been 
obtained using the background spacetime arena or the notions of space and time of the 
background metric. Such notions are unavailable at the singular point of the spacetime 
manifold. Unless the quantum theory of geometry or gravitation, whatever that means, is 
invoked the problem of accounting for the usual physical phenomena remains if any genuine 
spacetime singularity were to causally affect its exterior. In other words, the visible or naked 
spacetime singularity becomes a genuine problem for the classical theories of physics. 

It is well-known ||^ that solutions with visible or naked spacetime singularities can be 
obtained in the General Theory of Relativity as a held theory of gravitation. Construct¬ 
ing such a spacetime geometry and equating the corresponding Einstein tensor with the 
energy-momentum tensor always produces a solution of the held equations of general rel¬ 
ativity. However, it remains to show why such solutions cannot be considered physical in 
some dehnite sense of the term, ie. we require some principle on the basis of which visible 
spacetime singularities can be considered to be unphysical. This is the problem of the Cos¬ 
mic Censorship Hypothesis (CCH). Moreover, since the classical theories can be considered 
to be complete theories, in some appropriate sense, such a principle may be expected to be 
classical in nature or content. 

Thus one of the most important but as-yet unsolved problems of the classical theory of 
general relativity is undoubtedly that of the CCH. In fact, this problem is also important for 
the fundamental theories of physics since it directly concerns the nature of the gravitational 
held in the strong held limit. It is the same hypothesis that is at the helm of the very 
existence of Black Holes and hence Black Hole Physics, Astrophysics and Thermodynamics 
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as well [^. As a result it is certainly important to either prove or disprove the expectation 
that this principle is of a classical nature. 

Unfortunately, it has been extremely difficult to arrive at the precise and generally agree¬ 
able formulation of the CCH let alone prove it. This is primarily because of the conceptual 
nature of the problem itself. The statement of the CCH is that the spacetime singularity 
be not visible to any legitimate physical observer in a generic, realistic, physical situation 
such as gravitational collapse [Q. We emphasize that this is simply an expectation based 
on our inability to treat a genuine spacetime singularity in a classical theory. It necessarily 
follows that if any genuine spacetime singularity were to causally affect its exterior then the 
external observer would be unable to account for the physical phenomena in its vicinity. 
The classical predictability would then be in jeopardy if the singularity were visible to the 
external observer. 

As regards the statement of the hypothesis, the real issue then is about making the terms 
such as generic, physical or realistic precise so as to make some mathematically provable 
or disprovable statement for this principle of Cosmic Censorship. For example, we could 
demand that the collapsing matter satisfy some energy conditions/equation of state, that 
the spacetime not possess any special symmetry, that it should be stable to perturbations 
and so on. However, it has proven to be rather difficult to translate such demands into any 
generally agreeable mathematically precise statements. 

In the absence of any precise formulation of the CCH and hence its proof, we can, at 
best, look for examples of naked or visible singularities. From such examples, we can hope 
to sharpen the meaning of various terms in the statement of the CCH. 

In the recent past, many workers have followed this approach. However, we note 
that none of these examples can be considered to be a genuine counter-example to the CCH 
since all can be considered to be special in some or the other sense. Nonetheless, it is true 
that such examples can be used to provide some dehnite meaning to the statement of the 
censorship hypothesis. 

To illustrate the above point, we note that the locally naked singularity of the Vaidya-de 
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Sitter metric can be shown |]^ to be the same as that of the Vaidya metric. Then, we 
can argue that the asymptotic flatness of the solution does not manifest, in some sense, in 
the occurrence or not of the naked spacetime singularity, ie. the asymptotic observer has 
no role to play in the formulation of the censorship hypothesis. Therefore, this example 
unequivocally supports the formulation of the strong CCH due to Penrose |^. Essentially, 
the Strong CCH demands that the space-time singularity be not visible to any observer. 

In this paper, we considered the shear-free collapse of a spherical star with heat flux 
using the space-time of the metric ( p.l|) . Moreover, our model, based on the metric (p.l|) , can 
accommodate the realistic feature that the star can possess an extended atmosphere similar 
to what is observed with the real stars. Our results show that the singularity resulting in 
the stellar collapse modelled using the metric ( p.l|) is not visible to any observer if the initial 
density prohle of the star is non-singular similar to what is observed with real stars. The 
strong cosmic censorship hypothesis states that the singularity be not visible to any observer. 
Hence, our results are in complete accord with the Strong CCH for the non-singular density 
prohle of the stellar model considered here. 
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APPENDIX A: 


In this appendix, we discuss the energy conditions as obtainable for the considered space- 
time. (For an elegant analysis and further details of energy conditions applicable to imperfect 
fluids, see [pd]] .) 

We recall the dehnitions of the energy conditions for the sake of completeness: 

• Weak energy condition-. TabW^-W^ > 0 where is a timelike, future-directed 4-vector. 
It implies that the energy density as measured by any observer is positive. 

• Dominant energy condition-. Fa = TabW^ must be future-directed, timelike or null for 
any timelike, future-directed W°^. It implies that the speed of energy flow of matter 
must be less than the speed of light for every observer. 

• Strong energy condition-. 2 TahW^-W^ + T > 0 for any timelike, future-directed, unit 
4-vector W^- where T is the trace of Tab. 


We may note here that the Dominant energy condition implies the Weak energy condition, 
that the Strong energy condition can be violated in the presence of negative pressure, that 
is, when TabW^-W^ is negative and that the investigation of energy conditions is an algebraic 
problem that leads to a search of the roots of a polynomial of degree 4 in a four dimensional 
spacetime. 

From 1^, we also recall here the corollary for the shear-free or dynamic-viscosity-free, 
spherically symmetric spacetimes. We note that the metric ( 0 ) describes a shear-free 
spacetime. 

Fluid that undergoes shear-free motion or has vanishing coefficient of dynamic viscosity 
satishes 


• the Weak energy condition iff 


p + p>2q 

(Al) 

p — p 4- A > 0 

(A2) 
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the Dominant energy condition iff, in addition to (|A1|) , we have 


p — p>0 p — 3p + A>0 

• the Strong energy condition iff, in addition to ([^), we have 

2p + A > 0 


(A3) 


(A4) 


where A = \J{p + pY — 

In general, these conditions are to be imposed on the spacetime under consideration. 
However, it is easy to verify that these conditions are satisfied for the spacetime of the 


metric (p.l|) when 0 < a < 1 for the simple equation of state p = ap. 
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